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Abstract. The detection of the signature of dipole gravity modes has opened 
the path to study the solar inner radiative zone. Indeed, g modes should be the 
best probes to infer the properties of the solar nuclear core that represents more 
than half of the total mass of the Sun. Concerning the dynamics of the solar core, 
we can study how future observations of individual g modes could enhance our 
knowledge of the rotation profile of the deep radiative zone. Applying inversions 
on a set of real p-m.ode splittings coupled with either one or several g modes, we 
have checked the improvement of the inferred rotation profile when different error 
bars are considered for the g modes. Moreover, using a new methodology based 
on the analysis of the almost constant separation of the dipole gravity modes, 
we can introduce new constraints on solar models. For that purpose, we can 
compare g-mode predictions computed from several models including different 
physical inputs with the (?-mode asymptotic signature detected in GOLF data 
and calculate the correlation. This work shows the great consistency between the 
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signature of dipole gravity modes and our knowledge of p-modes: incompatibility 
of data with a present standard model including the Asplund composition. 



1. Introduction 

Solar modeling still needs improvements as constraints are not strong enough 
to correctly reconstruct for instance, the observed sound speed profile as well as 
the rotation profile in the solar core. 

With the release of the new abundances of Asplund et al. (2005), the discrep- 
ancy between the sound speed and density profiles of standard models including 
this new composition and the ones inferred from acoustic (p)-mode observations 
increased, especially in the radiative zone (Turck-Chieze et al. 2004a; Bahcall 
et al. 2005; Guzik et al. 2005; Mathur et al. 2007; Zaatri et al. 2007). Many 
attempts have been made to explain this discrepancy by attributing it to opaci- 
ties, equation of states, or by changing the composition in Ne or the microscopic 
diffusion rates (see Basu et al. 2008; Guzik et al. 2008 in these proceedings). 
As for the rotation profile, it is well established down to O.SRq (Couvidat et al. 
2003; Garcia et al. 2007) with larger and larger uncertainties if we go deeper 
into the solar core. Besides, the limitation of the p modes to bring information 
on the solar core rotation has been shown with inversions of data sets using real 
error bars (Garcia et al. 2008b). Indeed, by using more low-degree p modes 
(up to 3.4 ^Hz), we manage to reconstruct artificial rotational profiles down to 
0.2Rq (see Fig.[l]). 
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Figure 1. Results of the inversion of a simulated rotation profile (blue) using 
only p modes with three different sets of modes (see Garcia et al. 2008b). 

Gravity {g) modes would be the best data to improve our knowledge on the 
dynamics and the structure of the solar core. Today, we are almost 'touching' 
their individual detection and with a few more inches, we would be able to 
grab them. Here is a non exhaustive list of works dedicated to the search for g 
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modes: detection of one candidate at 220 /iHz (Garcia et al. 2001; Gabriel et 
al. 2002, Turck-Chieze et al. 2004b; Mathur et al. 2007), the component of this 
candidate followed in time with VIRGO (Garcia et al. 2008; Jimenez & Garcia 
2008 these proceedings), the detection of their asymptotic properties (Garcia et 
al. 2007). Moreover, new techniques are still being developed (Appourchaux 
et al. 2008; Garcia et al. 2008; Grec et al. 2008 these proceedings). In this 
work, we investigated what can be inferred about the solar core from g-mode 
detection. 



2. On the Dynamics of the Inner Radiative Zone 

To estimate what specific information on the core rotation is added by g modes, 
modes for which we have today candidates and that are likely to be unequivocally 
observed in the future, we used here an inversion code based on the Regularized 
Least Square methodology (Eff-Darwich &; Perez-Hernandez 1997). We created 
an artificial rotational profile, which is a step-like profile (see Figj2|. From this 
profile, we calculated artificial splittings and we inverted them after attributing 
them observational error bars. We have considered two sets of g modes: one 
with the current i^-mode candidate (D2 and D3) and another one with 8 low- 
degree g modes (£ = 1, n = -2 to -5 and £ = 2, n = -3 to -6) (D4 and D5), and 
for both cases uncertainties values of 75 (D2 and D4) and 7.5 nHz (D3 and D5). 
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Figure 2. Inversions for the simulated profile (blue). In the four sets of 
data, p modes used are as set 1 of Fig. [l] (see Mathur et al. 2008). 



Fig. [2] shows that with one g mode and a very large error bar, we do not have 
any information on the behavior of the rotation in the solar core. Then, either 
by increasing the number of modes or by decreasing the error bars, we start to 
have some indication of the actual simulated rotation rate in the core. However, 
such small amount of modes is not enough to reconstruct the profile. We only 
have one or two points below O.IR© in the inversion. We need more g modes 
to constrain the profile at these depths. For the region between 0.1 and 0.2 Rq, 
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low-degree high-frequency p modes are needed as their inner turning points are 
in this region. 
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Figure 3. Inversions of real data with only p modes and adding the g-mode 
candidate (from Mathur et al. 2008). 



Then, we apphed the inversion to real data using only p modes and also adding 
the g-mode candidate with the previous values of uncertainty. As different sce- 
narios were proposed to explain this candidate, two values of splitting have been 
used: 300 and 600 nHz. Whatever the value of the splittings is, the result is 
stable and we have the same behavior down to O.2R0. If we take into account 
the error bars, the rotation is well defined down to O.17R0. However, the rota- 
tion profile below this radius cannot be determined as long as we do not detect 
several g modes individually. 



3. On the Physical Inputs in Solar Models 

In this section, we tried to see how the detection of asymptotic properties could 
provide us some information on a few physical inputs in solar models. We know 
that g-mode periods are periodically spaced in the Power Spectrum Density 
(PSD). The Fourier transform of this PSD (PS2) has a large peak (APi for 
the modes i=l) corresponding to this periodicity (see for example Garcia et 
al. 2008c). Using the same methodology as the one developed in Garcia et 
al. (2007), we applied it to seven solar models with different physical inputs. 
Depending on the models used, ^-mode frequencies can vary up to 5;uHz. This 
information cannot be used as long as g modes are not detected individually. 
But the asymptotic value of APi is within a one-minute range depending on the 
physical input of the model (Mathur et al. 2007). 

The method consists in filtering only the region delimited by the peaks of the 
APi and its first harmonic in the PS2 (Fig. |4]left) so that we can reconstruct 
the signal corresponding to this region of the spectrum. This method is very 
powerful as it reconstructs the components of the modes (Fig. E] right) and it 
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is sensitive to the rotation. Indeed, in the PS2, we cannot distinguish between 
a model with a sohd rotation and a model with a core rotating faster than the 
rest of the radiative zone. Thus, as well as testing the dynamics of the core, we 
can check the sensitivity to some physical processes. Two parameters are used: 
the physical input such as the chemical composition, the microscopic diffusion 
or the diffusion in the tachocline, and the rotation profile which is a step-like 
profile with different values of rotation rate in the core. 




Period (hours) 

Figure 4. Reconstruction of signal: on the left (top), PS2 for 2 models with 
different rotation profile, as a function of period that is filtered on the peak 
of APi and its first harmonic (blue regions). On the left (bottom), signal 
reconstructed for both models. On the right, simulated signal (black peaks) 
and signal reconstructed for 2 rotation profiles (red and blue) (extracted from 
Garcia et al. 2007). 



Once we have reconstructed the signal corresponding to different solar models 
and to the GOLF data, we calculate the correlation rate between these models 
and the data. Fig. [5] shows the correlation rates for the seven models as a 
function of the rotation profile. For a given rotation profile, two groups can be 
distinguished: the first one with higher correlation rates for the models using 
the old abundances of Grevesse et al. (1993) and another group of models with 
the abundances of Asplund et al. (2005) with very low correlation rates (even 
negative) and the model without microscopies diffusion. Thus, GOLF data are 
more compatible with models S, and with the old abundances than the other 
models. Then, a rotation profile with a high rate from 3 to 5 times the rate in 
the rest of the radiative zone gives higher correlation rates. 

4. Conclusions 

On one hand, a theoretical work was carried to see how individual g-mode detec- 
tion could improve our knowledge on the rotation of the solar core. We showed 
that with one or a few g modes, we can have some information on the rotation 
rate in the core. But to better reconstruct the rotation profile in the core, we 
need more modes. On the other hand, we used the detection of asymptotic 
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Figure 5. Correlation rates between reconstructed signal of GOLF data 
and seven solar models as a function of the rotation profile (from Garcia et 
al. 2008c). 



properties to extract information on the structure of solar models. The method- 
ology used is sensitive to the structure for a few models and also to the rotation 
profile. Moreover, the results obtained are in agreement with p-mode results as 
GOLF data are more compatible with standard models using old abundances. 
We still have to keep in mind that these standard models are not complete and 
the discrepancy is not necessarily related to the abundances themselves. How- 
ever, to improve our knowledge on the solar burning core, we need to detect 
these g modes: individually, we would be able to constrain the dynamics and 
globally, we would have some constraints on physical inputs. 
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